The one-pion-exchange three-nucleon force and the A y puzzle 
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We consider a new three-nucleon force generated by the exchange of one pion in the presence 
of a 2N correlation. The underlying irreducible diagram has been recently suggested by the au- 
thors as a possible candidate to explain the puzzle of the vector analyzing powers A y and iTn for 
nucleon-deuteron scattering. Herein, we have calculated the elastic neutron-deuteron differential 
cross section, A y , iTn, T20, T21, and T22 below break-up threshold by accurately solving the Alt- 
Grassberger-Sandhas equations with realistic interactions. We have also studied how A y evolves 
below 30 MeV. The results indicate that this new 3NF diagram provides one possible additional 
contribution, with the correct spin-isospin structure, for the explanation of the origin of this puzzle. 
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The A y puzzle (or, more appropriately, the puzzle of 
the vector analyzing powers) is probably the most fa- 
mous of the open problems in 3iV scattering at low en- 
ergy. The problem with this observable has been ob- 
served quite early at the Tokio & Sendai Conference 
(Few-Body XI, 1986) B, since at that time the first 
reliable Faddeev calculations with realistic 2N potentials 
were becoming available, thanks in particular to the em- 
ployment of separable expansion methods which trans- 
form the 3 A'' scattering equations of Alt, Grassberger, 
and Sandhas (AGS) into an effective, multichannel 
Lippmann-Schwinger equation. This method of calcula- 
tion has been pushed forward to obtain accurate results 
particularly by the Graz group || . Since then, various al- 
ternative methods of solution of the 3A?" scattering equa- 
tion have been developed and tested 0] , and outstanding 
progresses have been made in the computational tech- 
niques in order to: 1) include three-nucleon forces (3NF) 
in the 3N scattering equations 2) treat explicitly 

the A dynamics in the 3N system M; 3) provide a com- 
bined description of the 3A^ dynamics with Coulomb, re- 
alistic 2N, and phenomenological 3N forces M. 

The puzzle was confirmed by these new approaches, 
and it turned out that the existing 2ir-3NF |p|-pd| pro- 



vided a too small effect for A y |4||5|,|8j , and not always in 
the right direction. In absence of new 3NF's that could 
explain the puzzle, and since the 3A^ A y is rather sensitive 
to the 3 Pj NN phase shifts, it was concluded in Ref. |l^] 
(see also references therein) that such phases and the as- 
sociated NN potentials derived from modern phase-shift 
analysis must be modified at low energy. These modifica- 
tions can be achieved without affecting appreciably the 
2N data because the low-energy 2A~ observables cannot 
resolve the 3 Pj phases uniquely due to the Fermi- Yang 
ambiguities. However, as has been argued in Ref. |jl3| , 
it is not possible to increase the 3A" A y with reasonable 
changes in the NN potential, hence additional 3NF's of 
new structure have to be considered. Recently, an at- 



tempt has been made ]14[ ] using a purely phenomenolog- 
ical 3NF of spin-orbit type, constructed ad hoc to affect 
only the triplet-odd states of 2N subsystem. 3NF terms 
of pion-range/short-range form |l5| have been reconsid- 
ered lately from the point of view of Chiral Perturbation 
Theory Jl6| (%PT), which predicts for these terms a non- 
negligible role. Qualitatively, it was found that these 
terms somewhat affect A y , but a quantitative conclusion 
could not be derived. 

It is our intention to indicate here a possible solution of 
this puzzle in terms of a new 3NF proposed recently by 
the authors jlTj . This force is generated by the one-pion- 
exchange diagram when one of the two nucleons involved 
in the exchange process rescatters with a third one while 
the pion is "in flight" . The underlying diagram has been 
derived starting from a formalism |18|] devoted to the 
explicit treatment of the pion dynamics in the 3N system. 
The resulting dynamical equation resembles a Faddeev- 
AGS equation, but entails in its inner structure the full 
complexity of the underlying four-body (ttNNN) system. 
A gradual procedure to project out the pion degrees of 
freedom has been discussed in Ref. where it has 

been shown that this formalism leads to irreducible 3NF 
diagrams, and includes in particular the 3NF we will 
evaluate here. 

3NF diagrams of the type derived in Ref. [[17| and ana- 
lyzed here are not new in the literature pffl-E2|, but they 
have been discarded in modern few-nucleon calculations 
because of the presence of a cancellation effect which has 
been observed in Refs. and discussed later from 

the point of view of effective chiral Lagrangians |25| . This 
cancellation involves meson retardation effects of the it- 
erated Born term, and the irreducible diagrams gener- 
ated by sub-summing all time orderings involving the 
combined exchange of two mesons amongst the three 
nucleons. However, this cancellation is incomplete (Tt]] 
and generates a three-nucleon force through a subtrac- 
tion term of the type H-v" with the 2N i-matrix being 
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pushed fully off-the energy shell because of the presence 
of the pion. Obviously, some ambiguities and model de- 
pendencies should be expected, since the subtraction in- 
volves the t-matrix in a region which is difficult to access 
and constrain, e.g., by 2N data. In addition, there are 
also possible model dependencies on how this "imper- 
fect" cancellation should be specified, since the input 2N 
potential itself could possibly include - maybe in some- 
what hidden way - meson-retardation correction effects. 
Ideally, this new 3NF term should be constructed con- 
sistently with the specific aspects and details of the given 
2N interaction. 

The explicit expression of this force (ItJ is 
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The full 3-/V-F is given by summing over the cyclic per- 



mutations of the nucleons, V 3N — V 3N 



vi N 



vi N . 



The momenta p,q represent respectively the Jacobi co- 
ordinates of the pair "12" , and spectator "3" , while E is 
the 3./V energy. We have set the pion-nucleon coupling 
constant to the traditional value f% NN / '(47r) — 0.078, and 
have employed standard form-factors of monopole type 
to describe the effective, composite nature of the meson- 
baryon coupling: 



firNN\Q) = firNN 



Al- 



K + Q 2 



(2) 



The chosen 2N model interaction implicitly determines 
the value of pion-nucleon cut-off in V 3N , e.g. = 1.7 
GeV for the Bonn B potential. The transferred momen- 
tum Q = q' q enters also in ui^ — yjm 2 + Q 2 . ty 
denotes the subtracted i-matrix between nucleons 1 and 
2, defined according to the prescription 



t 12 (p, p'; Z) = c{Z) t 12 (p, p'; Z) - v 12 (p, p') 



(3) 



Other details can be found in Rcf . |l7j . 

In addition, we have introduced here the effective pa- 
rameter c(Z), which is the only adjustable quantity of 
this 3NF. This parameter represents an overall correc- 
tion factor for the far-off-the-energy-shell 2N i-matrix 
entering this 3NF diagram. Ideally, c(Z) should be one 
for a 2NF model able to provide a reliable extrapolation 
of the t- matrix down to Z rj —160 MeV. However, none 
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FIG. 1. The A y puzzle, for nd scattering at 3 MeV (Lab). 
Calculations with the Bonn B potential (dots), and with the 
high-rank BBEST potential (lines). With the BBEST poten- 
tial there are two calculations, one obtained using the (non 
separable) 2-dim approach, the other with the separable 1-dim 
algorithm. The curves are not distinguishable. Data (trian- 
gles) from Ref. §§. 

of the existing 2N i-matrices can guarantee such extrap- 
olation since they are all constrained by experiments at 
the deuteron pole and at Z > 0. Furthermore, c(Z) 
might also correct for possible model dependencies on 
how this imperfect cancellation manifest itself, as already 
observed at the beginning of this communication. On 
general grounds, one expects that with increasing energy 
in nd scattering, the factor c{Z) should drift towards one, 
since the off-shell 2N t-matrix in V 3N approaches grad- 
ually the energy region with experimental constrains. 

To calculate the nd scattering observables below 
threshold, we have used the high-rank BBEST potential 
as two body input ||,^6| ■ With this separable representa- 
tions of the Bonn B potential, it is possible to solve accu- 
rately the Faddeev-AGS scattering equations, and obtain 
results comparable (with errors less than 1%) to those ob- 
tained from a direct solution of the Faddeev equations us- 
ing the original potentials as input. As an example, Fig. [j] 
shows the results we have obtained for A y at 3 MeV with 
the BBEST potential and with the original Bonn interac- 
tion. There are three curves since the calculations have 
been performed using both the separable (1-dimensional) 
algorithm and the non-separable (2-dimensional) method 
based on spline interpolation and Pade approximants, 
but the lines are practically indistinguishable. Similar 
tests have been made also before |27],^8| on various occa- 
sions. Triangles represent the nd experimental data from 
Ref. |§. 

The three-nucleon forces can be incorporated in the 
scattering equation in a relatively simple way if the 2N 
input potential is of finite rank. We sketch the proce- 
dure for a rank one case: once the separable 2N i-matrix 
is given, t = |<7i)t(<7i|, and the (anti)symmetrized AGS 
equation have been rewritten in the Lovelace form, 

X u = Z ll +Z ll rX ll , (4) 

the above 3NF is incorporated into this one-dimensional 
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FIG. 2. 3NF effects on the differential cross section and 
analyzing powers, for nd scattering at 3 MeV. The dashed 
lines are results with no 3NF. Solid line exhibits 3NF effects 
as discussed in the text. 

integral equation with the driving term calculated as fol- 
lows 



Zn = (gi\G P\ gi ) + {gi\G V^G \ 9l ) 



(5) 



The first contribution is the standard 2N driving term, 
with Go and P being the free Green's function and the 
cyclic/anti-cyclic permutator, respectively, while the sec- 
ond takes into account the effects of the 3 AT force. This 
procedure is consistent with the formalism developed in 
Ref. 19 to include irreducible pionic effects in the 3 A" 
dynamics, and at the same time it corresponds to the es- 
tablished perturbative procedure to include 3NF effects 
in the separable AGS equation fBOfl . 

The results obtained are shown in Fig. |[ for the differ- 
ential cross-section, and the analyzing powers Ay, iTn, 
T 2 o, T 21 , and T 22 . The panel exhibits the results obtained 
with the BBEST potential, where a rank-4 expansion has 
been used for all states with j < 2, aside for the states 
(rank 5), 3 S 1 - 3 D 1 (rank 6), and 3 P 2 - 3 F 2 (rank 5). 

The results without the 3A^F are shown by the dashed 
lines, and are practically indistinguishable from the cor- 
responding results for the Bonn B potential. The solid 
curve considers the additional contribution of the 3A" 
force expressed in Eq. ([!]) . We have used the subtraction 
method of Eq. (^), with t and v given by the Bonn B in- 
teraction while the parameter c(Z) has been set to 0.73 
for this energy. Results obtained with the PEST/Paris 
potential are similar. 

In Fig. ^ we show how the puzzle evolves above the 
break-up threshold, up to 30 MeV. In the top panel we 
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FIG. 3. Same as in previous figure for A y at 10, 18, and 
29.6 MeV. 

compare our theoretical calculations with the experimen- 
tal data at 10 MeV 31 1, obtaining the correct repro- 



duction for A y with the 3NF when c(Z) = 0.735. In 
the middle panel comparison is made with pd data at 18 
MeV |52|. Comparison of nd calculations with pd data 
is somewhat questionable because of the perturbations 
introduced by the Coulomb field. However it is known 
that, aside for the angles in forward direction, one of the 
main effects of the Coulomb field in A y can be approxi- 
mately reproduced by comparing the charged data with 
calculations performed at energies lowered of about 0.5- 
0.7 MeV. For this reason the data are compared with 
calculations at 17.3 MeV. (However, care must be exer- 
cized in interpreting this fact as a Coulomb slow-down 
effect [p3| ) At this energy we obtain the reproduction 
of the observable with c(Z) = 0.773. The botton panel 
shows the results obtained at 29.6 MeV, compared with 
the corresponding data at the same energy J34j (trian- 
gles), and with pd data at 30.2 MeV |35| (squares). The 
solid line has been obtained with c(Z) = 0.81. It is ev- 
ident that as the energy increases, the c(Z) parameter 
shifts slowly towards one, as expected. 

We checked also how the triton binding energy is af- 
fected, and found that here the effects are relatively 
small. With the highest possible rank, i.e., with a rank 
5 representation in all states with j < 2, except the cou- 
pled states 3 S'i- 3 i^i (rank 6), and 3 P 2 - 3 F 2 (rank 7), the 
BBEST+3NF result is -8.137 while the corresponding 
2NF result is -8.090 (MeV). 

We observe that some unavoidable approximations en- 
tered in this study. For instance, in obtaining Eq. (0) 
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we have ignored nucleonic recoil effects and have divided 
the subtracted t- matrix by the pion mass, instead of 
to get simpler expressions in partial waves. In addition, 
there might be a possible 3NF contribution of shorter 
range for the exchange of a p-meson; this term whould 
then counteract the one-pion-exchange 3NF, at least in 
the tensor part. Other uncertainties are related to the 
perturbative treatment of the pion dynamics in the AGS 
equation ]j~9fl . Finally, uncertainties about the fully off- 
the-energy-shell extrapolation of the NN i-matrix en- 
tering in this 3NF forced us to introduce the effective 
parameter c(Z). 

As discussed in Refs. |]l7| , |l9| , the 3NF contribution 
analyzed in this study is just one class of irreducible dia- 
grams generated by the pion dynamics, and in a more 
complete analysis also other classes of 3NF diagrams 
should be taken into account. Forces of the type TM ||, 
Brasil Jl(| , Urbana , belong to another class and they 
complement the 3NF whose effects have been calculated 
here separately. Some of these 3iVF"s give small correc- 
tions to A y , but not always in the right direction |p|,|l6|; 
their relevance, however, appears to be greater elsewhere 
(e.g., the binding energy of the triton). There are also 
additional 3NF terms of shorter range that might con- 
tribute, and in particular those obtained from x?T ]T^ ] 
appear very promising in providing an additional correc- 
tion to A y , possibly in the right direction. 

To summarize and conclude: we have evaluated here 
for the first time a new "pionic" effect in the 3N system. 
The effect is a natural consequence of a recently devel- 
oped theory [jl8],[l9| for the combined 7r-3iV dynamics in 
the 3N system, and has been recast into a 3NF term of 
new structure by the authors p"7| ] . The underlying 3 NF 
diagram complements the extensively discussed 2ir-3NF 
diagrams, and this complementarity shows up in the way 
this force affects the 3A^ observables: while the 2ir-3NF 
terms have a large contribution on the 37V binding energy 
and little effects (in the considered energy range) on the 
vector analyzing powers, we have shown here that this 
new 3NF term greatly modifies in particular these two 
spin observables, and has the potential to provide in full 
the solution of the A y puzzle. Conversely, we checked 
also that the same force produces smaller changes for 
the triton binding energy. Since both effects are clearly 
needed for describing the low-energy behaviour of the 3N 
system, it will be important to investigate at this point 
what will be the effect of the combined treatment of these 
two forces. 
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